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ABSTRACT 

We present a deuterium abundance analysis of the line of sight toward the white dwarf WD 2211—495 
observed with the Far Ultraviolet Spectroscopic Explorer (FUSE). Numerous interstellar lines are de¬ 
tected on the continuum of the stellar spectrum. A thorough analysis was performed through the simulta¬ 
neous fit of interstellar absorption lines detected in the four FUSE channels of multiple observations with 
different slits. We excluded all saturated lines in order to reduce possible systematic errors on the col¬ 
umn density measurements. We report the determination of the average interstellar D/O and D/N ratios 
along this line of sight at the 95% confidence level: D/O = 4.0 (±1.2) x 10 -2 ; D/N = 4.4 (±1.3) x 10 _1 . 
In conjunction with FUSE observations of other nearby sight lines, the results of this study will allow a 
deeper understanding of the present-day abundance of deuterium in the local interstellar medium and 
its evolution with time. 

Subject headings: ISM: abundances - ISM: clouds - cosmology: observations - ultraviolet: ISM - stars: 
individual (WD2211-495) - (stars:) white dwarfs 


1. INTRODUCTION 

Deuterium is believed to be produced in appreciable 
quantities only in primordial Big Bang nucleosynthesis 
(BBN) and destroyed in stellar interiors (e.g., Epstein 
et al. 1976); it is thus a key element in cosmology and 
in galactic chemical evolution (see e.g. Vangioni-Flam & 
Casse 1995; Prantzos 1996; Scully et al. 1997). The pri¬ 
mordial abundance of deuterium is one of the best probes 
of fish 2 , the baryonic density of the Universe divided by 
the critical density. The abundance of deuterium relative 
to hydrogen (D/H) is expected to decline during Galactic 
evolution at a rate that is a function of the star formation 
rate; standard models predict a factor of 2 to 3 decrease 
in the deuterium abundance in 15 Gyrs (see e.g., Galli 
et al. 1995; Tosi et al. 1998). Hence, any abundance of 
deuterium measured at any metallicity should provide a 
lower limit to the primordial deuterium abundance. This 
picture is essentially constrained by deuterium abundance 
measurements at look-back times of ~ 14 Gyrs (primordial 
intergalactic clouds), 4.5 Gyrs (protosolar), and 0.0 Gyrs 
(interstellar medium). Although the evolution of the deu¬ 
terium abundance seems to be qualitatively understood, 
measurements of D/H at similar redshifts show some dis¬ 


persion and indicate that additional processes may be im¬ 
portant in controlling the abundance of deuterium. That 
fact has led to the development of non-standard models, 
which propose, for example, larger astration factors (e.g., 
Vangioni-Flam et al. 1994) or non-primordial deuterium 
production (see e.g., Lemoine et al. (1999) for a review). 

Up to now, the interstellar medium is the astrophysical 
site that has allowed the most comprehensive investiga¬ 
tions of deuterium abundances. Deuterium has been ob¬ 
served in the interstellar medium using different methods: 
radio measurements of its 92 cm hyperfine transition (e.g., 
Blitz & Heiles 1987; Chengalur et al. 1997), observations 
of deuterated molecules (e.g., Lubowich et al. 2000; Ferlet 
et al. 2000), Balmer series analyses (Hebrard et al. 2000), 
and Lyman series absorption (see Moos et al. 2001). Of 
these, the most accurate measurements have been ob¬ 
tained through Lyman absorption-line observations in the 
far-ultraviolet (far-UV) spectral range. By observing hy¬ 
drogen and deuterium directly in their atomic form, far- 
UV Lyman series absorption-line measurements provide 
accurate column density determinations that are not de¬ 
pendent on ionization or chemical fractionation effects. 
The first measurement of the abundance ratio (D/H)ism 
was reported by Rogerson & York (1973) for the line of 
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sight to (3 Cen, using Copernicus'. (D/H)ism= 1.4 ± 0.2 x 
10 -5 . Since then, many other (D/H)ism measurements 
have been performed using different instruments [Coper¬ 
nicus, IUE, IMAPS, HST) for other sight lines, and the 
values obtained show significant dispersion around the 
above value. 

For example, an average value (D/H)ism= 
1.50(±0.10) x 10 -5 (la) has been derived for the Lo¬ 
cal Interstellar Cloud (Lallement & Bertin 1992) by Lin- 
sky (1998) from the comparison of 12 nearby sight lines, 
but studies of several lines of sight revealed values outside 
this range ( e.g ., Laurent et al. 1979; York 1983; Allen et 
al. 1992; Vidal-Madjar et al. 1998; Hebrard et al. 1999; 
Jenkins et al. 1999; Sonneborn et al. 2000). This dis¬ 
persion may result from spatial variations due to some 
unknown physical processes or underestimation of sys¬ 
tematic errors. There is still considerable debate over 
these two interpretations, and the final resolution of the 
issue may have implications for understanding the physics 
of the interstellar medium, as well as the baryonic den¬ 
sity inferred from D/H measurements. Measurements of 
D/H in the intergalactic medium are also sparse and do 
not agree with each other [e.g., Webb et al. 1997; Buries 
& Tytler 1998a; 1998b), so an accurate determination 
of the primordial deuterium abundance has also proven 
to be a difficult quantity to measure. Moreover, recent 
studies of the anisotropy of the cosmic microwave back¬ 
ground (CMB), which permits evaluations of the baryonic 
density independent of those obtained through deuterium 
measurements (see, e.g., de Bernardis et al. 2000; Jaffe et 
al. 2001), imply higher values of fish 2 than those implied 
by the abundance studies, corresponding to a primordial 
D/H value of the order of the (D/H)ism values (see how¬ 
ever de Bernardis et al. 2001; Stompor et al. 2001; Pryke 
et al. 2001). 

An accurate determination of the interstellar deuterium 
abundance is one of the main objectives of the Far Ul¬ 
traviolet Spectroscopic Explorer [FUSE), which was suc¬ 
cessfully launched on 1999 June 24 (Moos et al. 2000). 
In this paper, we present new measurements of deuterium 
abundances obtained with FUSE toward the white dwarf 
WD 2211—495. This paper is part of a series dedicated 
to FUSE measurements of deuterium interstellar abun¬ 
dances toward BD+28° 4211 (Sonneborn et al. 2001), 
WD 1634-573 (Wood et al. 2001), WD 0621-376 (Lehner 
et al. 2001), G191-B2B (Lemoine et al. 2001), HZ43A 
(Kruk et al. 2001), and Feige 110 (Friedman et al. 2001). 
A major goal of this program is to determine to what ex¬ 
tent and on what scales variations in the D/H ratio occur. 
Moos et al. (2001) present an overview of these first results. 

White dwarfs are ideal targets for the measurement 
of the interstellar deuterium abundance (Lemoine et 
al. 1996); they may be chosen close to the Sun (so that 
the column densities are not too high and the velocity 
structure of the line of sight is not too complex), and they 
exhibit relatively smooth UV continua (which allow detec¬ 
tions of lines from many species). Prior to this series of pa¬ 
pers, measurements have been published for (D/H)ism to¬ 
ward three white dwarfs: G191—B2B (Lemoine et al. 1996; 
Vidal-Madjar et al. 1998; Sahu et al. 1999), Sirius-B 
(Hebrard et al. 1999), and Feige 24 (Vennes et al. 2000). 


The FUSE observations of WD 2211—495 and the data 
processing are presented in Sect. 2, and the details of the 
analysis in Sect. 3. The results are reported in Sect. 4 and 
discussed in Sect. 5. 

2. OBSERVATIONS AND DATA PROCESSING 

WD 2211-495 (RE J2214-491) was identified by the 
ROSAT wide field camera all-sky survey as an extreme ul¬ 
traviolet source and was classified as a new white dwarf 
star (Pounds et al. 1993). Holberg et al. (1993) analyzed 
the optical and high-resolution International Ultraviolet 
Explorer [IUE) spectra of WD 2211—495 and concluded 
it was a hot DA white dwarf showing traces of heavy el¬ 
ements such as C, N, O, Si, and Fe. Later, Holberg et 
al. (1994) and Werner & Dreizler (1994) discovered the 
presence of Ni in the co-added IUE spectrum, confirm¬ 
ing that WD 2211—495 was one of the most metal-riclr 
white dwarfs known. Barstow et al. (1998) determined the 
atmospheric parameters of the star by considering non- 
LTE metal-line-blanketed model atmospheres and found 
T eff = 62,000 K and log# = 7.2. WD 2211-495 has a 
mass M ~ 0.53 M 0 according to optical spectroscopic 
analyses performed by Marsh et al. (1997a), Vennes et 
al. (1997), and Finley et al. (1997). Its photometric dis¬ 
tance is d ~ 53pc (Holberg et al. 1998). Table 1 summa¬ 
rizes relevant sight line and atmospheric parameters for 
WD 2211-495. 

2.1. Observations 

WD 2211—495 was observed in 1999 and 2000 as part of 
the FUSE photometric calibration program and the FUSE 
Science Team D/H program. Seventeen observations con¬ 
sisting of 188 individual exposures for a total exposure 
time of 93 800 sec (~ 26 hours) were used in the present 
study. All the observations were obtained in histogram 
mode, through one of the three slits available (Moos et 
al. 2000): the large (LWRS), medium (MDRS), or high 
resolution aperture (HIRS). Table 2 shows the log of the 
observations used in this analysis. 

2.2. Data processing 

The one-dimensional spectra were extracted from the 
two-dimensional detector images and calibrated using 
slightly different versions of the CALFUSE pipeline (see Ta¬ 
ble 2), depending of the version available in the archives 
for a given observation. For a given slit (LWRS, MDRS, 
or HIRS), data from each channel and segment (SiCIA, 
SiC2B, LiFlA, LiF2B, etc.) were co-added separately, af¬ 
ter wavelength shift corrections of the individual calibrated 
exposures. Wavelength shifts between individual spectra 
were typically a few pixels for exposures of a given obser¬ 
vation and a few tens of pixels between exposures of dif¬ 
ferent observations. Because a given wavelength falls on 
different pixels in different exposures, co-addition of the 
numerous exposures results in a considerable reduction in 
fixed-pattern noise (Sahnow et al. 2000) by acting like a 
random FP-Split procedure (Kruk et al. 2001). 

Exposures with strong airglow emission, large spectral 
resolution variations, or large wavelength stretching 3 were 
not included in the sums. The last two criteria were ap¬ 
plied in particular to the two observations performed in 
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1999, early in the mission. In order to properly co-add 
different exposures for a given segment, the same wave¬ 
length calibration file (wave*009) was used in CALFUSE for 
all of the observations. That led us to re-run CALFUSE in 
the version 1.8.7 (which includes wave*009) for the obser¬ 
vations which were archived in an older version which did 
not include this calibration file (see Table 2). The spectral 
resolution in the final spectra ranges between ~ 13000 and 
~ 18500, depending on detector segment and wavelength 
(see Section 4.1). 

All of the spectra were binned to three pixel samples 
(the line spread function (LSF) is about 10 pixels wide). 
The final LWRS spectra are plotted in Figure 1; similar 
datasets were obtained for the MDRS and HIRS slits. A 
sample of a spectrum (SiClB with LWRS) is shown in 
Figure 2 with the D I, O I, and N I lines identified. 

3. DATA ANALYSIS 

3.1. Method, overview 

For a given transition of a given element, each inter¬ 
stellar cloud along the sight line produces an absorption 
line that can be modeled with a Voigt profile. In addition 
to atomic parameters [taken from Morton (1991; 1999)] 
a Voigt profile is defined by four cloud parameters: the 
radial velocity v of the cloud (in kms -1 ), the column 
density IV) of the element i (in cm' 2 ), the temperature 
T of the gas (in K), and its micro-turbulent velocity 
£ (in kms - ). FUSE does not have sufficient spectral 
resolution to resolve the individual velocity components 
along the WD 2211—495 sight line. This implies that 
FUSE cannot resolve the fine shape of individual veloc¬ 
ity components and that it cannot reveal the presence of 
closely spaced multiple components. We do not have high- 
resolution ground-based interstellar data at our disposal 
for this star. 

The analysis has been done using the profile fitting pro¬ 
cedure Owens.f, developed by one of us (M.L.) and the 
FUSE French Team. A subset of the data has been inde¬ 
pendently fitted using a second procedure (see Sect. 4.2.2), 
developed by another of us (B.E.W.). The reader should 
refer to Wood et al. (2001) for a full description of this 
fitting procedure, and we will concentrate on describing in 
detail the results obtained with Owens . f. 

Owens. f models Voigt profiles using % 2 minimization. 
One of the characteristics of this procedure is that it is 
able to fit simultaneously several lines in different spectral 
windows, a spectral window being a part of a full spectrum 
(see Figures 3-6). Thus, this procedure finds the best solu¬ 
tion compatible with numerous spectral lines based on the 
assumption that all lines considered yield the same values 
for v, T, and £ for a given component, and that all lines for 
a specific component of element i give the same value for 
Ni. The width of the lines combines both parameters T 
and £, which can only be separately determined if several 
elements with different masses are simultaneously fit 4 . 

Data at the edge of each detector segment are inaccu¬ 
rate (see Figure 1) and were not used in the fit. In the 
end, 32 different absorption lines (D I, O I, N I, Fe II, 
Si II, and P II) were modeled, which were observed in dif¬ 


ferent segments and spectra, resulting in a total of 115 
simultaneously fitted lines. Remaining ions like Ar I and 
N II were modeled separately in extra fits. No H 2 lines 
were detected. The 115 fitted lines are located in 65 dif¬ 
ferent spectral windows, extracted from the different spec¬ 
tra, obtained through the HIRS, MDRS, and LWRS slits, 
on channels SiCl, SiC2, LiFl, and LiF2. Plots of the 65 
spectral windows are shown in Figures 3-6. Atomic pa¬ 
rameters of the D I, O I, and N I lines used are reported 
in Table 3. 

3.2. Free parameters 

All the cloud parameters (v, T, £ and all the JVj’s) were 
free to vary in the fit. We also allowed the continua shape, 
the spectral shift and the LSF width to vary from one 
spectral window to the next. These parameters are a pri¬ 
ori unknown or poorly known, and this uncertainty is likely 
to introduce systematic effects in the final result. We take 
into account this uncertainty by treating these unknowns 
as free parameters in the fit, so that their uncertainty is 
included in the y 2 variations (see Sect. 3.6). 

The continua were fit by 3 rd to 7 th order polynomials 
(see Figures 3-6), depending on the spectral region. All 
the parameters of the polynomials were free. The spectral 
shift of each window 5 (i.e., its wavelength zero point) was 
free to allow for possible inaccuracy of the FUSE wave¬ 
length calibration. 

All Voigt profiles were convolved with a simple Gaussian 
LSF. The spectral resolution of FUSE is not accurately de¬ 
termined, but it is known to depend on the segment and 
wavelength. We let the width of the LSF vary freely for 
each spectral window. Note that fits with a double Gaus¬ 
sian LSF were also performed (see Sect. 4.2.2). 

3.3. Unsaturated lines 

Since we are primarily interested in abundances, we fo¬ 
cused on the determination of column densities. We chose 
to use only unsaturated lines in the fits to reduce uncer¬ 
tainties resulting from saturation effects and other system¬ 
atic errors. Saturated lines are not very sensitive to the 
column density. Since the lines which we use are on the 
linear part of the curve of growth, we can measure the to¬ 
tal column densities for the line of sight even though the 
individual clouds are not resolved. All the H I interstellar 
lines detectable in the FUSE spectral range are located 
on the flat part of the curve of growth for this sight line; 
hence, no accurate H I column density could be deter¬ 
mined. We thus determined D/O and D/N ratios instead 
of D/H. In order to fit the Lyman D I lines without fitting 
H I lines with Voigt profiles, the blue wings of the H I lines 
were fit with polynomials (see e.g., Figure 5, left panel) to 
provide the continuum for D I absorption (see, however, 
Sect. 4.2.2). 

3.4. Stellar model 

We compared the FUSE spectrum of WD 2211—495 to 
a stellar model in order to verify that the interstellar 
lines were not blended with stellar lines. To this end, 


4 The 6-value, which is the width of a line, is given by b~ = 2AL _|_ where k is the Boltzmann constant and m the mass of the element. 

5 Actually, the spectral shift was free for all windows but one in order to avoid degeneracy in a parameter space with the radial velocity of the 
absorber. 
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we computed a non-LTE metal line-blanketed model at¬ 
mosphere using the programs TLUSTY/SYNSPEC (Hubeny 
& Lanz 1995). We used the atmospheric parameters de¬ 
termined by Barstow et al. (1998) (see Table 1). We per¬ 
formed an abundance analysis using the FUSE and IUE 
data and then computed a synthetic spectrum incorpo¬ 
rating these abundances (Chayer et al. 2001). The com¬ 
parison between the stellar model and the FUSE spectrum 
allowed us to exclude the O I A950.88A and Si II A989.87A 
interstellar lines that are blended with the P IV A950.66A 
and N III A989.80A stellar lines. In agreement with the 
DA-type classification of the WD 2211—495, no He II ab¬ 
sorption is seen in the FUSE spectrum ( e.g ., A1084.9A or 
A992.3A), and hence such absorption should be absent at 
the intervening H I Lyman lines. 

3.5. Zero flux level 

To properly deduce abundances from interstellar absorp¬ 
tion lines, it is extremely important to precisely evaluate 
the zero flux level. Although scattered light in FUSE is 
low (Moos et al. 2000; Sahnow et al. 2000), some residual 
flux is visible on the bottom of the H I Lyman lines, as one 
can see in Figures 1 and 2. This residual flux (at a level of 
few percent of the continuum) is probably due to the scat¬ 
ter of adjacent flux from the continuum due to the wings 
of the LSF. This interpretation agrees with the fact that 
this small amount of flux is not detected at wavelengths 
lower than the Lyman break, where there is no flux for 
scatter. Since the residual flux is so low, for each spectral 
window we tuned the zero flux level to the signal at the 
core of the nearest Lyman line, where it is readily measur¬ 
able because these lines are totally saturated. The diffuse 
light is a few percent of the continuum, and its variations 
as a function of wavelength seem to be of the same order 
of magnitude, or even lower. Fits with a double Gaussian 
LSF were also performed (see Sect. 4.2.2). 

3.6. y 2 scaling and Ay 2 method 

The % 2 of the final fit was 2424.2. The difference be¬ 
tween the number of independent spectral bins used in the 
fit and the number of free parameters (cloud, polynomials, 
shifts, and LSF parameters) was 1764. For this number of 
degrees of freedom, the reduced y 2 is 1.37. This value is 
larger than 1.00, probably because of a small underestima¬ 
tion of the tabulated errors for each spectral bin. These 
errors are computed by the extraction pipeline which does 
not take into account all possible instrumental artifacts. 

For each estimated parameter, e.g. N(D I), error bars 
were computed from an analysis of the y 2 variation. We 
pegged that parameter at various trial values, and for each 
trial value we ran an extra fit, allowing all the other pa¬ 
rameters to vary freely. We used these minima to define 
the confidence interval. Scanning the parameter in that 
way for different values (up to at least Ay 2 = 50), we ob¬ 
tained the value of y 2 as a function of this parameter, and 
we derived the 1, 2, 3, 4, 5, 6, and 7 a error bars using the 
standard Ay 2 method. The error bars at these different 
confidence levels have been taken into account to produce 
the final 2cr uncertainties reported here. The Ay 2 curves 
as a function of column density are plotted in Figure 7. 

However, since our best fit y 2 value lies far above what 
would be expected for the corresponding number of de¬ 


grees of freedom, and since that change is due to inaccu¬ 
rate estimates of the noise array, we rescaled our y 2 and 
the corresponding Ay 2 values in the following way. At 
the 95% confidence level, the smallest value that one can 
draw from a y 2 distribution with 1764 degrees of freedom 
is 1667.5, and consequently we increase the noise variance 
on each pixel by a factor 2424.2/1667.5 ~ 1.45 (2424.2 cor¬ 
responds to our best fit value), or equivalently we divide 
our y 2 and Ay 2 by a factor 1.45. Note that this is the 
most conservative choice at the 95% confidence level for 
the rescaling parameter. For example, for the 5 a contour 
in y 2 one should now look for differences Ay 2 = 25 x 1.45 
using our original error array, which corresponds to the 
standard Ay 2 = 25 for the correctly rescaled error array. 
This rescaling, similar to that presented by Lemoine et 
al. (2001), has increased in a conservative way the error 
bars derived below. 

4. RESULTS 

4.1. Parameters values 

Our column densities measurements are reported in Ta¬ 
ble 4. We derived the following abundance ratios: 

D/O = 4.0 (±1.2) x 10 -2 

D/N = 4.4 (±1.3) x 10 -1 

with 2 er error bars. These error bars were computed from 
the individual errors on the column densities listed in Ta¬ 
ble 4, assuming that they are uncorrelated. Although the 
noise may not follow a Gaussian distribution, we checked 
on smooth parts of the spectrum that the tabulated error 
bars scaled as described above (Sect. 3.6) effectively in¬ 
clude more than 95% of the individual points of the spec¬ 
trum. Therefore, the quoted “2 cr” error bars can be con¬ 
sidered as a 95% confidence level, statistic and systematic 
effects being taking into account (see Sect. 4.2). 

A formal temperature and turbulent velocity measured 
through the fit are T = 3400 ± 2500K and £ = 3.0 ± 
1.0km s -1 . These values should be taken with caution be¬ 
cause several velocity components may be present and the 
spectral lines are not resolved. It is also a combination 
between turbulent velocity and velocity dispersion of the 
possible different clouds. In addition, these values were 
obtained from fits of unsaturated lines, which are not very 
sensitive to temperature and turbulent velocity. 

For a given channel, the relative shifts between differ¬ 
ent spectral windows were found to be relatively constant, 
with a dispersion below 10 kms -1 and a standard devi¬ 
ation of ~ 3 kms -1 . Shifts between spectral windows of 
different channels were larger (up to 130 kms -1 ). 

The full widths at half maximum obtained were 9.9±1.4, 
10.9 ± 1.3 and 10.9 ± 1.2 pixels respectively for spectral 
windows from HIRS, MDRS, and LWRS (values are av¬ 
erages and standard deviations on all the corresponding 
windows). As expected, the spectral resolution is slightly 
higher for HIRS than for MDRS and LWRS. All the widths 
found were in the range 8-13 pixels. Note that for 14 
spectral windows it was not possible to determine the LSF 
width because of the low signal-to-noise ratio; in these 
cases, the width was tuned to that obtained in the adja¬ 
cent spectral windows. 
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4.2. Systematic effects 
4.2.1. Overview 

Error bars on column densities computed with the Ay 2 
method (Sect. 3.6) include statistical and many possible 
systematic effects. The statistical aspect is obvious since 
the statistical error on each pixel is taken into account in 
the x 2 computation. 

Several approaches were used in order to reduce system¬ 
atic errors. Some effects were treated as free parameters 
(Sect. 3.2). Some systematic effects are likely to be de¬ 
creased by fitting simultaneously several different lines for 
a given element. It is possible that some of these effects 
affect our fits but the use of numerous different lines for a 
given element (see, for example, Table 3 for D I, O I, and 
N I) averages out any possible bias. In the same spirit, ef¬ 
fects which are totally or partially linked to instrumental 
effects are likely to be decreased by having data from mul¬ 
tiple segments, multiple observations, and multiple slits 
fitted simultaneously. It is unlikely that one of these ef¬ 
fects affect in the same manner a given line observed with 
what are essentially different instruments (different seg¬ 
ments are used, through different slits, on different parts 
of the detectors). Finally, many possible systematic er¬ 
rors in the column density measurements are reduced by 
analyzing only unsaturated lines. 

We list here the possible systematic effects which are 
treated using these four approaches: 

• Free parameters: 

— continuum; 

— spectral shifts; 

— line spread function width. 

• Simultaneous fit of several lines for a given element: 

— line blending; 

— oscillator strengths uncertainty; 

— continuum; 

— zero flux level determination; 

— fixed-pattern noise effects; 

— line spread function effects; 

— flux calibration; 

— wavelength calibration. 

• Simultaneous fit of several segments: 

— continuum; 

— zero flux level determination; 

— fixed-pattern noise effects; 

— line spread function effects; 

— flux calibration; 

— wavelength calibration. 

• Unsaturated lines: 

— number of interstellar components; 

— temperature of the clouds; 

— turbulence of the clouds; 

— shape and width of the line spread function. 


4.2.2. Discussion of the effects and tests 

In this section we discuss all the effects we investigated 
and we present some tests that we carried out in order to 
check the reliability of the error bars. 


• Line blending 

Line blending with unknown interstellar, stellar, or air- 
glow features may affect the fit. However, such features 
are not present systematically for all lines of a given el¬ 
ement. One exception is the He II stellar lines, falling 
systematically near the interstellar D I lines. However, no 
He II absorption is seen in the WD 2211—495 spectrum in 
the unblended lines at 1084.9A or 992.3A. 


• Oscillator strengths uncertainty 

Oscillator strengths (/) of D I lines are well known. How¬ 
ever, this is not the case for some lines of other elements, 
such as O I or N I. If some /-values were found to be in¬ 
accurate, the fit of that line would stand out as a poor fit 
and the line would be rejected ( e.g ., O I 1026.47A). By 
fitting each line individually, we checked that each tabu¬ 
lated /-value is consistent with the others (none of them 
present an abnormal tabulated /-value). The range in 
column densities from fits to the individual lines retained 
in the sample presented a ~ 20% dispersion around the 
value obtained with all the lines. Maybe some oscillator 
strengths are slightly inaccurate, but these errors should 
be small and average out in the final result. 


• Continuum placement 

The continua of the interstellar lines is given by the stellar 
continuum, possibly modified by instrumental effects (for 
example, one can see in Figure 1 the bump near 1165A on 
LiFl, which is due to an instrumental artifact). For each 
spectral window, the continuum fitted by a polynomial is 
near the actual continuum but is probably not a perfect 
model. However, when the continuum is fitted indepen¬ 
dently for each line, in spectra obtained in each segment 
and each spectrograph aperture, it is unlikely that the av¬ 
erage placement will be systematically higher or lower than 
the actual continuum. 

We tested the possible systematic effect of the degree 
n of the polynomials by performing fits with values of n 
ranging from 1 to 14. Column densities obtained were con¬ 
sistent within the 2 a error bars. Thus, there is no strong 
dependence of the fit on the degree of the polynomials used 
to fit the continua. 

Around the D I lines, the continuum is affected by the 
blue wing of the H I Lyman line. In order to check the D I 
continuum levels found using polynomial approximations 
to the wings of the H I Lyman lines, we performed addi¬ 
tional fits with Voigt profiles used to describe the H I lines. 
The value of IV(D I) remained within the 2cr error bar, so 
we infer that continuum levels at the D I lines are not 
strongly dependent upon how the H I lines are modeled. 
We note that the H I column density is not constrained 
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well by these fits since the H I Lyman lines lie on the flat 
part of the curve of growth. 

Finally, we divided the spectra by the stellar model we 
computed (Sect. 3.4); fits of these normalized spectra pro¬ 
duce column densities included within the 2 a error bars. 
Thus, no significant systematic effect linked to the contin¬ 
uum was found. 

• Zero flux level determination 

The zero flux level may be a function of wavelength. Si¬ 
multaneous fits of different lines of a given element average 
out possible inaccuracies in the zero flux level determined 
from the bottom of the Lyman lines (Sect. 3.5). 

We made extra fits in two extreme cases: doubling that 
correction (increasing the zero flux level by comparison 
with the best fit) and letting it remain zero. We also made 
fits where the zero flux level was a free parameter on each 
spectral window, with all windows fitted simultaneously. 
All these fits produced column density values again within 
the 2 a error bars. 

• Fixed-pattern noise effects 

Fixed-pattern noise is able to create false features that 
may artificially increase or decrease the intensity of an ab¬ 
sorption line. The comparison of results for any single line, 
measured in different segments and apertures, allowed us 
to find lines affected by fixed-pattern noise. In the final 
fit, no lines were found to be in strong disagreement with 
the other ones so these effects are small and average out 
in the fit. 

• Datasets 

In addition to fitting all the data simultaneously as de¬ 
scribed above, we also made fits to partial subsets of the 
total dataset presented in Table 2. These subsets were 
chosen to find possible time-dependent effects, segment to 
segment, channel to channel or detector to detector dif¬ 
ferences, and slit (HIRS, MDRS, LWRS) effects such as 
optical transmission, fixed-pattern noise, and airglow con¬ 
tamination (airglow intensity depends of the size of the 
aperture). All of the column densities found from these 
partial fits remained within the 2 a error ranges, so no 
significant effects were found. As expected, error bars ob¬ 
tained with only one of these subsets of data were larger 
than the ones obtained with the complete dataset. 

• Reduction procedure 

Fits were made from observations reduced with different 
spectra extraction, different wavelength calibration files, or 
different extraction windows on the two-dimensional im¬ 
ages of the detectors. We found no significant effects due 
to the calibration uncertainties. 

• Selection of lines 

In order to identify possible effects on a given line (such 
as fixed-pattern noise right at the wavelength of a line 
or an erroneous /-value), we also fit each line indepen¬ 
dently of all others. No abnormal results were found for 
any line when compared to the composite fit results. In 


addition, this allowed us to identify saturated lines to be 
excluded from the final fit. For example, fits to the O I 
lines AA948.69, 971.74 and 1039.23 (each line taken inde¬ 
pendently), with line strengths /A larger than ~ 5A, had 
large uncertainties and overestimated the value of N(0 I) 
derived from weaker transitions. We rejected the saturated 
lines using this objective criterion. 

• Selection of the species 

The D I, O I, N I, Fe II, and Si II lines were fitted as¬ 
suming that all these species were in the same interstellar 
(neutral) cloud. The same radial velocity, turbulence, and 
temperature were assumed for all these species. This as¬ 
sumption may be incorrect if the species do not co-exist in 
the same regions along the sight line, which might be true 
if a significant amount of ionized gas is present. In order 
to test the possible systematic effect due to this assump¬ 
tion, we made extra fits including all these species but one. 
We also performed a fit including all the species but Fe II 
and Si II, which are the most likely ions to trace ionized 
gas. For all these fits, results for each column densities 
remained within the 1 a error bars. 

• Line spread function 

We used a simple Gaussian LSF in the fitting procedure 
and adopted a small zero flux level shift to account for 
weak, residual light contributed probably by a broad com¬ 
ponent of the LSF (Sect. 3.5). Note that the residual light 
does not affect the weak lines used in this study in an ap¬ 
preciable way. The robustness of this approach was tested 
on a simulated spectra made with a double Gaussian LSF; 
they were fitted with a simple Gaussian LSF and a tuned 
zero flux level, and the result was satisfactory (see below). 

We also attempted to fit all spectral windows with a 
unique LSF width. Column densities obtained remained 
within the 2 a error bars for LSF full widths at half max¬ 
imum included in the range 9 12.5 pixels; Ay 2 of the 

9 and 12.5 pixel LSF solutions were large (> 200). Note 
that the minimum y 2 was obtained for a 10.8 pixel width 
LSF. However, this y 2 was greater by ~ 200 than that 
obtained with the LSFs left free to vary from one spectral 
window to the next, showing that as expected, the LSF 
does depend on the wavelength, the segment, and the slit. 

Finally, we performed fits with double Gaussian LSFs. 
We made tests using the LSF established by Wood et 
al. (2001), and other tests using free double Gaussian LSFs 
(Lemoine et al. 2001). Fits obtained with double Gaus¬ 
sian LSFs gave the same column densities within the 2er 
error bars. 

• Temperature and turbulent velocity 

The values obtained for T and / are inaccurate (Sect. 4.1). 
We note, however, that the the impact of the uncertainty 
on the column densities is weak; for temperatures and tur¬ 
bulent velocities within 0 - 20000K and 1.5 - 7.0km s -1 , 
the column densities are consistent within the 2 a error 
bars. Hence, no effects on column densities due to T or £ 
are significant, as expected from the fact that only unsat¬ 
urated lines are fit. 
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• Number of interstellar components on the sight line 

Fits were made assuming only one interstellar component 
on the line of sight, even if it is likely that several compo¬ 
nents are present along this ~ 50 pc sight line. Since we 
use only unsaturated lines, assuming only one interstellar 
component should have no effect on the measurement of 
the total column densities. Fits obtained with up to five 
interstellar components along the line of sight gave the 
same total column densities within the 1 a error bars. 

• Simulation 

In order to test the various approaches of analyzing the 
interstellar lines from the FUSE spectra, an artificial spec¬ 
trum was created for the FUSE Team members (Moos et 
al. 2001). This simulated spectrum includes two inter¬ 
stellar clouds, He II lines in the stellar continuum, random 
noise, fixed-pattern noise, high frequency ringing, dithered 
wavelengths, background, and a double Gaussian LSF de¬ 
pending on the segment. This simulated spectra was fit in 
a blind manner using the method presented above (includ¬ 
ing only one interstellar component and simple Gaussian 
LSF). All the estimated column densities were in agree¬ 
ment with the inputs in the simulation within the 2 a error 
bars. This validates our approach and led us to a better 
understanding of the errors and the possible systematic 
uncertainties that might be present in our analysis. 

• Method 

Finally, a subset of the data presented here has been inde¬ 
pendently fitted using the procedure presented by Wood et 
al. (2001). Both approaches gave similar results, in agree¬ 
ment within 2 a. The final values listed in Table 4 reflect 
the combined effort of these analyses. 

5. DISCUSSION 

The O I and Ar I column densities presented here 
agree with the first FUSE results reported by Jenkins et 
al. (2000). Agreement is poor for N I, for which the value 
14.02 ± 0.15 was reported by these authors. This early re¬ 
sult was obtained with a lower number of lines and a lower 
signal-to-noise ratio. In the current study, the deficiency 
of N I along this sight line with respect to the value in B 
stars appears not to be as large as reported by Jenkins et 
al. (2000). 

A value D/H ~ 1.5 x 10 -5 is obtained from the deu¬ 
terium column density reported in Table 4 and the H I col¬ 
umn density obtained from EUVE spectra: log./V(H I) = 
18.76 (Wolff et al. 1998). No error bar is available for 
this model-dependent H I estimate. We note that the far- 
UV spectrum of WD 2211—495 cannot be interpreted by 
using homogeneous H+He models (Marsh et al. 1997b), 
and that Holberg et al. (1998) quoted a different result, 
log./V(H I) = 18.3, based on fitting the Lyman a line ob¬ 
served with by IUE. 

Another way to estimate D/H is to use our present 
D/O measurement together with the accurate abundance 
of interstellar oxygen obtained by Meyer et al. (1998): 
O/H = 3.19(±0.14) x 1CT 4 . The low dispersion of O/H 
confirms that O I is a good tracer of H I in the Galac¬ 
tic disk. The Meyer et al. (1998) measurement was ob¬ 


tained from GHRS observations of 13 lines of sight rang¬ 
ing from 130 to 1500 pc, with most of them closer than 
500 pc. This study examines interstellar clouds that 
are more distant than WD 2211—495. However, we note 
that the O/H value given by Meyer et al. (1998) agrees 
with that measured in the more local ISM (Moos et 
al. 2001). Using the Meyer et al. (1998) O/H value to¬ 
gether with our new FUSE D/O measurement, we obtain 
D/H = 1.3(±0.4) x 10 -5 (2 a). Note that using the revis¬ 
ited value O/H = 3.43(±0.15) x 10 -4 from Meyer (2001), 
which adopt an updated O I oscillator strength, we obtain 
D/H = 1.4(±0.4) x 10" 5 (2 a). 

In the same spirit, Meyer et al. (1997) determined from 
GHRS observations toward 7 stars the interstellar ratio 
N/H = 7.5(±0.4) x 10 -5 . Adopting this value, our D/N 
value gives D/H = 3.3(±1.2) x 10 -5 (2a). We can prob¬ 
ably understand the difference between the O-based and 
N-based determinations of D/H in terms of the likelihood 
that N could be somewhat more ionized than H, whereas 
O and H are more strongly coupled to each other by 
charge exchange reactions. For example, Vidal-Madjar et 
al. (1998) have found toward G191—B2B that the O/H is 
perfectly consistent with the average of Meyer et al. (1998), 
while N/H is found to be a factor 2 smaller than the Meyer 
et al. (1997) ratio. In the case of WD 2211—495, this pos¬ 
sibility is supported by the fact that Ar I is well below 
its expected column density in relation to both N I and 
O I, using the abundance ratios of these elements in B- 
type stars as a reference (Cunha & Lambert 1992, 1994; 
Holmgren et al. 1990; Keenan et al. 1990). According to 
arguments presented by Sofia & Jenkins (1998), this con¬ 
dition indicates that a substantial fraction of the material 
is partially ionized by EUV photons that are allowed to 
penetrate the region. At the same time, if n e ?rni so 
that ordinary recombinations with free electrons dominate 
over charge exchange reactions of N II with H I, we might 
expect N I to behave in a manner similar to Ar I, i.e., 
n Nii/ftNi > nHii/nnu since the neutral forms of both 
N and Ar have photoionization cross sections larger than 
that of H. This inequality with the H ionization fraction 
is less likely to arise with O, since the charge exchange 
cross section of O II with H I (Field & Steigman 1971) is 
considerably larger than for N II with H I (Butler & Dal- 
garno 1979). An approximate indication of the expected 
magnitudes of the deficiencies of N I and Ar I is shown in 
Figure 2 of Jenkins et al. (2000). 

In any case, an accurate H I column density value is 
needed to obtain a reliable (D/H)ism measurement along 
this line of sight. Eventual observation of the Lyman a 
toward WD 2211—495 should provide this value since this 
line is not on the flat part of the curve of growth and there¬ 
fore the damping wings can be used to measure 7V(H I) 
more accurately. As of now, D/O and D/N estimates ap¬ 
pear to be more robust than D/H for the WD 2211—495 
sight line. 

Since H I, O I, and D I have nearly the same ionization 
potential, D/O ratio is an important proxy for the D/H 
ratio and its putative variations. D/O is very sensitive 
to astration, both from D I destruction and O I produc¬ 
tion (Timmes et al. 1997). If FUSE measurements of the 
D/O ratio toward other stars shows a relatively uniform 
value with low dispersion, it would argue in favor of D/H 
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and 0/H homogeneity in the interstellar medium. Indeed, 
the only other possibility would be that D/H and O/H 
vary precisely in the same way as to cause D/O to remain 
constant. It seems improbable since (i) O/H appears to 
be uniform in the interstellar medium over paths of sev¬ 
eral hundred parsecs, and (ii) astration should lead to an 
anti-correlation of D I and O I abundances. An upcoming 
survey of D I and O I absorption along about 10 local lines 
of sight observed by FUSE should help determine whether 
the D/O ratio varies significantly from one sight line to the 
others (Hebrard et al. 2001). Initial indications are that 
the ratio does not vary substantially (Moos et al. 2001). 

6. CONCLUSION 

We have presented FUSE observations of the 905-1187 A 
spectrum of the white dwarf WD 2211—495. The column 
densities of several ions were measured through simulta¬ 
neous fits to the numerous unsaturated absorption lines 
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Fig. 1. The eight segments of the FUSE WD 2211 —495 spectra obtained through the large aperture (LWRS). Each channel (SiCl, SiC2, 
LiFl, and LiF2) is divided in two segments, separated by a gap. Similiar spectra were obtained for the two other slits (MDRS and HIRS); in 
total, 8 x 3 = 24 co-added spectra were used for the analysis. The bump in the LiFl continuum at 1150A is an instrumental artifact. 



Fig. 2.— FUSE SiClB spectrum of WD 2211—495 observed through the LWRS slit. Positions of D I, O I, and N I interstellar absorption 
lines are indicated. 
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Fig. 5. — Same as Fig. 3, for MDRS LiFl (upper) and LiF2 (lower) spectral windows. 


Table 1 

Target Summary for WD 2211-495 


Quantity 

Value 

Reference 

Spectral Type 

DA 

1 

l 

345.79° 

2 

b 

-52.62° 

2 

d (pc) 

53 

3 

V 

11.71 

4 

U-B 

-1.21 

4 

B-V 

-0.34 

4 

Teff (K) 

62,000 

5 

log g (cm s -2 ) 

7.2 

5 

log(He/H) 

< -4.0 

1 

Trad (km s- 1 ) 

+33.9 

6 


References. — (1) Holberg et 
al. (1993); (2) Hog et al. (1998); (3) 
Vennes et al. (1997); (4) Marsh et 
al. (1997a); (5) Barstowet al. (1998); (6) 
Holberg et al. (1998). 
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Fig. 7.— Ax 2 curves for D I, O I and N I. Each point on these curves corresponds to an individual fit, made with all the parameters free, 
but the D I, O I, or N I column density fixed. Ax 2 = 25 corresponds to a 5 a error bar on the column density. Ax 2 = 36 corresponds to 6 a. 
A\' 2 = 49 corresponds to 7 a. x 2 are rescaled on these plots (see Sect. 3.6). 


Table 2 

Log of the observations. 


Obs. date 

Obs. reference 

'Y' a 

AT b 

1 v exp 

Aperture 6 

CALFUSE d 

1999.10.25 

M1030303 

4.2 

8 

LWRS 

1.8.7 

1999.10.31 

M1030304 

3.0 

6 

LWRS 

1.8.7 

2000.06.03 

M1030305 

5.2 

11 

LWRS 

1.7.3 

2000.06.29 

M1030306 

4.2 

7 

LWRS 

1.7.7 

2000.08.17 

M1030307 

5.3 

11 

LWRS 

1.7.7 

2000.10.24 

M1030308 

4.1 

9 

LWRS 

1.7.7 

2000.10.24 

M1030309 

5.1 

10 

LWRS 

1.7.7 

2000.10.24 

M1030310 

5.8 

8 

LWRS 

1.7.7 

2000.10.25 

M1030311 

6.1 

13 

LWRS 

1.7.7 

2000.10.25 

M1030312 

5.5 

10 

LWRS 

1.7.7 

2000.06.03 

P1043801 

16.5 

35 

MDRS 

1.8.7 

2000.06.02 

M1030201 

4.8 

10 

MDRS 

1.8.7 

2000.08.18 

M1030202 

4.8 

10 

MDRS 

1.7.7 

2000.10.24 

M1030203 

4.8 

10 

MDRS 

1.7.7 

2000.06.02 

M1030101 

4.3 

9 

HIRS 

1.8.7 

2000.08.18 

M1030102 

5.8 

12 

HIRS 

1.8.7 

2000.10.24 

M1030103 

4.3 

9 

HIRS 

1.8.7 

Total: 6 

17 obs. 

93.8 

188 




a Total exposure time of the observation (in 10 3 s). 

b Number of individual exposures during the observation. 

C LWRS, MDRS, and HIRS are respectively large, medium and narrow 
slits. 

d Version of the pipeline used for spectral extraction. 
e Totals for all datasets used in this study. 
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Table 3 


D I, O I AND N I LINES INCLUDED IN THE FIT. 


A a 

/ b 

# c 

A a 

/ b 

# c 

A a 

f h 

# c 



D I 




01 



N I 


1025.4434 

7.91 

x 10" 2 

9 

930.2566 

5.37 x 10 -4 

5 

965.0413 

4.02 x 10“ 3 

5 

972.2723 

2.90 

x 10" 2 

6 

929.5168 d 

2.36 x 10" 3 

5 

964.6256 

9.43 x 10“ 3 

6 

949.4847 

1.39 

x 10" 2 

3 

925.4420 

3.54 x 10" 4 

3 

963.9903 

1.48 x 10 -2 

3 

937.5484 

7.80 

x 10" 3 

1 

924.9520 d 

1.59 x 10" 3 

6 

954.1042 

6.76 x 10 -3 

6 

930.4951 

4.82 

x 10" 3 

5 

921.8570 d 

1.19 x 10~ 3 

3 

953.9699 

3.48 x 10“ 2 

6 

925.9737 

3.18 

x 10" 3 

3 

919.9080 

1.78 x 10" 4 

3 

953.6549 

2.50 x 10“ 2 

6 

916.9311 

7.23 

x 10 -4 

2 

919.6580 d 

9.47 x 10" 4 

4 

953.4152 

1.32 x 10 -2 

6 





916.8150 d 

4.74 x 10" 4 

2 

952.5227 

6.00 x 10 -4 

2 








952.4148 

1.70 x 10 -3 

2 








952.3034 

1.87 x 10 -3 

2 



total: 

29 


total: 

31 


total: 

44 


a Wavelength in vacuum at rest (in A). 
b Oscillator strength. 

c Number of independent lines included in the fit (observed on different segments and through 
different slits). 

d Triplet structure used (see Morton 1991; 1999) 


Table 4 


Total interstellar column densities. 


species 

logAl(cm 2 ) a 

species 

logA7(cm 2 ) a 

D I 

13.94 (±0.10) 

Si II 

14.00 (±0.15) 

01 

15.34 (±0.08) 

Ar I 

12.82 (±0.10) 

Nl 

14.30 (±0.06) 

P II 

12.05 (±0.20) 

Fe II 

13.25 (±0.15) 

Nil 

14.6 (±0.3) 


2er error bars. 



